Stochastic process leading to wave equations in dimensions higher than one 
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Stochastic processes are proposed whose master equations coincide with classical wave, telegraph, 
and Klein-Gordon equations. Similar to predecessors based on the Goldstein-Kac telegraph process, 
the model describes the motion of particles with constant speed and transitions between discreet 
allowed velocity directions. A new ingredient is that transitions into a given velocity state depend on 
spatial derivatives of other states populations, rather than on populations themselves. This feature 
requires the sacrifice of the single-particle character of the model, but allows to imitate the Huygens' 
principle and to recover wave equations in arbitrary dimensions. 
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I. INTRODUCTION 

Can a wave, which is in general an essentially dy- 
namical process, be mapped onto a kinematic stochastic 
model like random walks? This question has a long his- 
tory and was approached from many different perspec- 
tives, ranging from pragmatic (numerical simulation of 
wave-related phenomena) to fundamental (interpretation 
of quantum mechanics) with many important incentives 
in between, such as effects of inertia in heat transfer, 
light propagation in turbid media, turbulence diffusion, 
etc. While random walks are natural underlying pro- 
cesses for parabolic equations of diffusion type, the con- 
nections between stochastic motion and hyperbolic wave 
equations are less obvious and often restricted to one spa- 
tial dimension (ID). Perhaps the best-known example 
is the telegraph equation f tt + \ ft = c 2 A/, which de- 
scribes propagation of waves in media with losses (with 
characteristic dissipation time r). In ID the equation 
can be readily derived from persistent random walk with 
constant speed and Poissonian velocity reversals (often 
referred to as the Goldstein-Kac process) (H-Q) but the 
same walk extended to higher dimensions does not evolve 
according to the telegraph equation 0, 0-0] ■ This fea- 
ture is generic and inherited in many related problems, 
in particular of mapping relativistic quantum wave equa- 
tions onto classical random walks [9l4lll| . Such mappings 
are typically designed in ID, and the extension to higher 
spatial dimensions requires the formal replacement of the 
space-variable derivative by the gradient, d/dx — > V. It 
was noted by many authors that this approach may be 
inconsistent since in general it is impossible to construct 
the random walk with desirable properties (governed by 
a master equation of desirable form) as a mere superpo- 
sition of independent one-dimensional walks. 

As will be discussed below, the difficulty of extending 
wave-particle isomorphism beyond ID is not related to 
stochastic nature of the random walk models, but rather 
originates from the inability for a single-particle motion, 



neither stochastic nor deterministic, to imitate the Huy- 
gens' principle in dimensions higher than one [l2j . On the 
other hand, there are often no reasons to restrict oneself 
to single-particle models, which cannot reproduce physi- 
cally relevant negative solutions anyway. In this paper we 
discuss a model a la Goldstein-Kac based on an ensemble 
of classical particles moving with a fixed speed and sub- 
jected to transitions between discreet allowed directions 
of motion. It is shown that transitions can be chosen in a 
form which allows to recover wave-like equations for the 
ensemble's distribution function in any dimension. 



II. WAVE EQUATION 

As well known, the classical wave equation in ID 
ftt = c 2 fxx can be mapped onto a totally determin- 
istic kinematic model. Thanks to the factorization 
(d/dt — cd/dx)(d/dt + cd/dx)f — 0, the general solu- 
tion can be written as the superposition of the functions 
f^(x,t), which satisfy the equations 



df+ 
dt 



df+ 
dx 



df_ 

dt 



dx 



(1) 



and can be interpreted as the distribution functions for 
independent particles, or for an ensemble of single par- 
ticles, moving freely with constant speed c in positive 
and negative directions, respectively. Eqs. ([1]) also 
can be written in terms the total distribution function 
f = f + + f~ and the current J = c (/+ — f~), 



dt 



dJ_ 

dx ' 



dJ_ 
~dt 



d£ 

dx 



(2) 



One might suggest that for higher dimensions the proper 
generalization of Eqs. J2j should read 



df -, dJ 

dt ~ V J ' dt 



-c 2 V/, 



(3) 



which indeed immediately gives the multidimensional 
wave equation for /, 
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d^£ 
dt 2 



c 2 A/, 



(4) 
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as well as the conservation law for the current vorticity 



0. 



(5) 



However, it is easy to see that the equations ([3]) do not 
follow from multidimensional generalization of Eqs. ([T]) 
and therefore cannot describe a single-particle motion in 
D > 1. 

It is instructive to illustrate the last statement explic- 
itly for two spatial dimensions (2D), assuming that par- 
ticles move with a constant speed c and can be in one of 
four velocity states (x, ±), (y, ±), corresponding to mo- 
tions along positive (+) or negative (— ) directions of the 
two Cartesian axes. Let f x (x, y, t) and f*(x, y, t) be the 
corresponding distribution functions. (From here on, we 
use subscripts a = x,y,z only to refer to vector com- 
ponents, not to partial derivatives). For freely moving 
particles there are no transitions between the states, so 
the two-dimensional version of Eqs. (JJ reads as follows: 



dt 



dfl 
dx 



9fl 
dt 



Tc- 



9Jl 
dy ' 



(0) 



Let f x = /+ + f x , f y = /+ + f y , and J x (x,y,t) = 
c (.fx - fx) and Jy{ x >V)t) = c (fy ~ f y ) are Cartesian 
components of the current J = iJ x + jJ y . Adding and 
subtracting Eqs. ([5]) one obtains 



df x _ d,J x 
dt dx 



dJ, 



c 2 dfx 

dt dx 



for particles moving along x-axis, and 



dfy 
dt 



dJy 

" dy ' 



dJy 

dt 



dfy 
dy ' 



(7) 



(8) 



for particles moving along y-axis. These equations do 
not form the closed system (j3|) for the total distribution 
f(x, y, t) = fx + f y and the current J, 



dj_ 
dt 
dJ_ 
~dt 



= -V-J, 



(9) 



:2)7 a f* , l d fv V ^ -c z Vf. 



dx dy 



Therefore the two-dimensional wave equation is not re- 
covered. 

Let us now generalize Eqs. ^ by allowing transi- 
tions between (x, +) and (x, — ) states with transition 
rate g x (x, y, t), and between (y, +) and (y, — ) states with 
transition rate g y {x, y, t), 



dfl 
dt 



dfl 
dx 



T9x, 



dt 



Tc^T 9y . (10) 



(Note that no transitions are still allowed between x and 
y states.) The presence of transitions does not violate the 
conservation of the number of particles, so the continuity 



equation df /dt — —V ■ J still holds, while the equation 
for the current now reads 



dJ_ 
~dt 



; - ,'dfr 2 \ , (dfy 2 



dx c ® x ) ^ V dy 



(hi 



(11) 



By choosing transition rates in the form 



9x 



c dfx 
2 dy 



the right-hand side of Eq. (fTTj) is completed to the gradi- 
ent, dJ /dt = — c 2 V/. Thus the system ([3]) is recovered, 
and therefore the total distribution / is governed by the 
two-dimensional wave equation Q. 

With transition rates (fT"2|). the partial motions along 
x and y axes become statistically coupled even though 
there are no direct transitions (x, ±) 4-> (y,±). Note 
that an equation for the particle distribution in a given 
state, say (x, +), 



dfl 
dt 



.dfl 
dx 



C dfy 

2 dx ' 



(13) 



has a form of the conservation law df x /dt = —dJ x /dx, 
where the flux </+ = c (/+ + hfy) is determined not only 
by the density of the given state /+ but also by local 
particles in other velocity states fy and /~. This is 
reminiscent of the Huy gens' principle, according to which 
every point of a wave front propagating with a speed 
c is the source of secondary waves that spread out in 
all directions with the same speed c. The resemblance 
is achieved at the expense of losing the single-particle 
status of the original ID model. It is easy to see that 
with transition rates (|12[) . proportional to derivatives of 
state populations, the partial distributions f x and f^ are 
not positively defined. Then Eqs. ([TO]) can not describe a 
single-particle motion, but should be interpreted as the 
equations for corresponding perturbations in an ensemble 
of particles. 

The generalization of the scheme for 3D is straightfor- 
ward and reads 



dfl 
dt 



df* 



a = x,y, z 



with 



c f dfy df z 
2 \ dx dx 



(14) 



(15) 



and similar expressions for g y and g z . The equations (|14[) 
and (fT5j) lead to the system ([3]) and therefore to the wave 
equation in 3D. 

A trick of coupling of partial motions with transition 
rates in the form (ITOI) or (fTSj) is quite generic and can be 
applied to derive other multidimensional wave-like equa- 
tions. Two more examples are presented below. 
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III. TELEGRAPH EQUATION 

A hybrid of the wave and diffusion equations, the tele- 
graph equation has a form 



dt 2 T dt 



c 2 Af 



(16) 



and describes the propagation of waves in media with 
losses. It is also often used as an approximation to treat 
diffusion processes beyond the overdamped limit and in 
a turbulent medium [l], [2, 8]. A one-dimensional version 
of Eq. fTo) can be derived as a master equation for the 
Goldstein-Kac stochastic process - a dichotomous persis- 
tent random walk in which a particle moves in II? with 
the velocity fluctuating between c and — c. Using the 
same notations as in the previous section, let f + (x,t) 
and f~(x,t) be the probability density for the particle 
moving to the right and to the left, respectively. Rever- 
sals of velocity are Poisson distributed and occurring with 
the rate l/2r. The processes is described by equations 



df+ 
dt 

at 



df+ 
dx 

df~ 
dx 



--if 
2t u 



2r 



-/"), 
-/")• 



(17) 



In terms of the total distribution function f = f + + f~ 
and the current J — c (f + — / _ ), the equations take the 
form 

df dJ dJ 1 T 2 df 

This leads immediately to the ID telegraph equation for 
/ 



dt 2 



1 



df 
at 



,0V 

dx 2 ' 



(19) 



and also for J and each components / + and / — . 

The isomorphism between persistent random walk and 
dissipative wave propagation, though very attractive 
from many points of view, does not go beyond II?. It is 
true that the multi-dimensional telegraph equation f| 16[) 
follows readily from the generalization of (CE3 



df 
dt 



~TTT = — V • J , 



dJ 1 ? 
dt t 



-c 2 V/. 



(20) 



However, these generalized equations do not follow 
merely from multi-dimensional extension of ID random 
walk (|17[) . For instance, such an extension for 21? has a 
form 



dfl 
dt 

df£ 
dt 

dt 
dt 



+c 



+c 



dfl 
dx 

df£ 
dx 

dy 

dj^_ 
dy 



^ (-3/+ + /- + fy), 



(21) 



-3f-+f++fy), 



•3/ a + fy 



fx), 

+ ^ ("3/„- + /+ + fx), 



1 



where f x = /+ + f x and f y = f+ + f y . Adding and 
subtracting lead to the following equations for the partial 
densities 



df* 
dt 

dfy 
dt 



dJ x 
dx 
dJ, 



-(.-fx + fy), 
T 



(22) 



and for the current components 



dJ x 
dt 
dJ„ 



dx 
2 df y 



Jx 1 

T 



dt dy 
Thus, for the total distribution function f = f x 
the current J = iJ x + j J y one obtains 

dj_ 
dt 



-irr = -V-J, 



(23) 



- f y and 



(24) 



dJ_ 

dt 



1 



-J = -c 2 



dfx , -/ df. 



dx 



+ 3 



dy 



? -c 2 Vf. 



This differs from Eqs. (f20j) . and therefore the two- 
dimensional telegraph equation is not recovered. 

Let us modify Eqs. (|21l) in precisely the same way as in 
the previous section. Namely, in addition to transitions 
with isotropic rate l/2r, let us introduce anisotropic 
transitions (x, +) <->■ (x, —) and (y, +) <-> (y, — ) with the 
rates g x and g y , respectively, 



dfl 

dt 
dfx 

dt 

dfl 
dt 

dj^ 
dt 



+c 



dfl 
dx 

df£ 
dx 

dfl 
dy 

dj^_ 

dy 



9.< 



g x 



f x 
ft 



9v + ^-(-m + fv 



2t 
1 



7»),(25) 

fx) J 



+ 9y + ^(-3fy-+f+ + fx) 



with g a given by Eq. ([12 



9x 2 



c dfy 
dx ' 



9y = g 



(26) 



c % 

In this case Eqs. (j2"2")l for df a /dt do not change, while 
Eqs.flM]) for dJ a /dt take the form 



dJ x 
dt 

dJ y 
dt 



dy 



(27) 



Thus both equations (|20|) are satisfied, from where the 
two-dimensional telegraph equation (fT6|) for / follows im- 
mediately. 

The generalization of the system (|25[) for 31? is 
straightforward, i.e the equation for f^ takes the form 

'' // = Tc^TQx+^i-Sft + fx* + fy + fz), 



2r 



dt 1 dx ' ^ ' 2r 

with f„ = /+ + /- and g x given by Eq. (JTS}. 
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IV. KLEIN-GORDON EQUATION 



For instance, the equations for take the form 



The combination of a substitution / = exp(— t/2r)ip 
and analytic continuation of the transition rate to the 
imaginary value l/2r — > — i A converts the telegraph 
equation (flB) into the Klein-Gordon equation 



c 2 AiP - A 2 ip. 



(28) 



This observation was exploited in many works discussing 
the possibility of the connection between random walks 
and quantum mechanics. Yet the concept of a complex 
transition rate does not look particularly attractive. 

An alternative stochastic model leading to the Klein- 
Gordon equation in ID without recourse to analytic con- 
tinuation was discussed in (Tlj . In this model transitions 
between two velocity states are governed by a "guiding 
field" E(x, t) which is coupled to the total particle distri- 
bution f(x,t) via the Poisson equation. Using the same 
method as in the previous sections, the model can be 
extended for 2D as follows: 



dfl 
dt 

dt 

dt 

dj^ 
dt 



= +c 



dx 
df£_ 
dx 

d Jl 
dy 

211 

dy 



9x 
9x 
9y 

9y 



1 



(29) 



v 



— aE„ 
2 y 



Here the transition rates g a are given by Eq. ([12 

_ C dfy 



c dfx 
2 dy ' 



(30) 



the field E = iE x + jE y satisfies the Poisson equation 



W-E = bf, 



(31) 



a and b are constants. For the total distribution / = 
fx + f y = fx + fx + fy + fy and the current J = 
c(/+ - f~)i + c(/+ - fy)j one obtains from Eqs. ([29j) 



dt ~ V J ' 



(32) 



dJ_ 
dt 



= -c 2 V.f + acE. 



Together with (jBTj) . these equations give the Klein- 
Gordon equation ([28]) for / with A 2 = a be. The possi- 
bility to interpret this approach in the spirit of the Dirac 
sea model was discussed in 



The generalization of this scheme for 3D is straight- 
forward and involves transition rates g a in the form (1151) . 



dfl 
dt 

df£ 
dt 

which gives 



= +c 



dfl 
dx 

df£ 
dx 



dfy 

dx 

dfy 

dx 



d_h 

dx 

d_h 

dx 



+ -aE x , (33) 



1 



- -^aE x 



dfx 
dt 

dJ x 
dt 



dJ x 
dx 

dx 



(34) 



acE x 



These and similar equations for y and z components lead 
to (|32p and therefore to the Klein-Gordon equations in 
3D. Note that if E is a conservative vector field, V x 
E = 0, then it follows from (I3"2l that the vorticity of the 
current V x J ( "spin" ) is a constant of motion. 



V. CONCLUSION 

In this paper we constructed systems of partial differ- 
ential equations of first order which can be interpreted 
in terms of stochastic motion of an ensemble of parti- 
cles and lead to wave equations for the particle distri- 
bution in any dimensions. There is nothing unusual, of 
course, in wave-like dynamics in a system of interacting 
particles. However, the presented model attempts to pic- 
ture waves not as a dynamical process but kinematically, 
assuming that the particles move with speed equal to 
the phase speed of the wave and represent medium ex- 
citations rather than medium constituents. Technically, 
the approach is based on a simple trick of "completing 
the gradient". Namely, an appropriate ID single-particle 
process is extended to higher dimensions supplemented 
with additional terms g a , which make the equation for 
the local current to involve the gradient of the total dis- 
tribution (e.g., dJ /dt = —c 2 Vf for the wave equation). 
The additional terms g a can be interpreted as transition 
rates between states with opposite velocity directions and 
make the model consistent with the Huygens' principle. 
The discreteness of velocity directions, assumed in this 
approach, may look physically artificial, yet technically 
it is essential. The model can be readily extended for 
any finite number of allowed directions, but not for a 
continuum of directions. 
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